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Modified theories of gravity have recentiy been studied by several authors as possibly viable alter- 
natives to the cosmological concordance model. Such theories attempt to explain the accelerating 
expansion of the universe by changing the theory of gravity, instead of introducing dark energy. 
In particular, a class of models based on higher order curvature invariants, so-called f(R) gravity 
models, has drawn attention. In this letter we show that within this framework, the expansion 
history of the universe does not uniquely determine the form of the gravitational action and it 
can be radically different from the standard Einstein-Hilbert action. We demonstrate that for any 
barotropic fluid, there always exists a class of f{R) models that will have exactly the same expansion 
history as that arising from the Einstein-Hilbert action. We explicitly show how one can extend the 
Einstein-Hilbert action by constructing a f(R) theory that is equivalent on the classical level. Due 
to the classical equivalence between f(R) theories and Einstein-Hilbert gravity with an extra scalar 
field, one can also hence construct equivalent scalar-tensor theories with standard expansion. 
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INTRODUCTION 

One of the great surprises in modern cosmology has 
been the accelerating expansion of the universe. The 
first sign of this unexpected behaviour came from dis- 
tant supernova observations (SNIa) in 1998, when it was 
discovered that far away supernovae are dimmer than ex- 
pected 0. The SNIa observations have since improved, 
notably by a number very distant supernova observations 
by the Hubble Space Telescope, and the main conclusion 
still stands: a critical density matter dominated universe 
is not consistent with the data. Such a universe, also 
known as the Einstein-deSitter (EdS) model, was also 
disfavoured by independent observations based on mea- 
suring the small fluctuations in the cosmic microwave 
background (CMB) @. Combining these two observa- 
tions along with a third independent observation based 
on the clustering of galaxies on very large scales 0], the 
present day cosmological model emerges. The cosmo- 
logical concordance model is a geometrically flat ACDM 
model, which is a critical density universe dominated by 
cold dark matter and dark energy in the form of a cos- 
mological constant. Furthermore, cross-correlating the 
observations of galaxy clusters with CMB measurements, 
one has seen evidence of the decay of gravitational poten- 
tials on very large scales Q, as predicted in the ACDM 
model. Although independent measurements are statisti- 
cally not very significant, combining measurements leads 
to a picture that is very well consistent with the concor- 
dance ACDM model Q. 

A natural alternative to adding a mysterious new form 
of energy to our physical picture of the universe, is to 
consider modifications of general relativity. This is also 
motivated by the fact that we only have precision ob- 



servations of gravity from sub-millimeter scales to solar 
system scales, which is very far from the present Hubble 
radius that is the scale relevant to the question of dark 
energy. 

A particular class of models that has drawn a signif- 
icant amount attention recentlyare the so-called f(R) 
gravity models (see e.g. 0, M H E3 

and references 

therein). They form a class of higher derivative gravity 
theories that include higher order curvature invariants 
in the gravitational action. Such an extension of the 
Einstein-Hilbert (EH) Lagrangian can be viewed natu- 
ral as there is no a priori reason why the gravitational 
action should be linear in the Ricci scalar R. Further- 
more, higher order terms can naturally appear in low en- 
ergy effective Lagrangians of quantum gravity and string 
theory. 

After f(R) theories were proposed as an alternative so- 
lution to the dark energy problem ^3|"E3> ^ was quickly 
realized that two major obstacles exist. Firstly, some of 
these theories exhibit instabilities that will render space- 
time unstable Q|. A second serious but more practical 
hindrance has been the difficulty in actually integrating 
the equations of motion numerically and comparing to 
observational data. 

Another twist in the tale is the question of frames and 
conformal transformations. While the Einstein and Jor- 
dan frames are mathematically equivalent on the classical 
level, the physical equivalence of the two frames at the 
perturbation and quantum level has been discussed in 
reference to f(R) theories. Possible restrictions arise due 
to non-standard gravitational effects constrained by fifth 
force experiments (see 0] and references therein). 

Recently, an inverse approach to the problem of inte- 
grating the complicated equations of motion arising from 
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a general f(R) action has been proposed in ^^|. Instead 
of specifying a particular form of f(R) one considers the 
inverse problem of reconstructing f(R) from the expan- 
sion history of the universe. 

In this letter we are also concerned with the inverse 
problem of reconstructing f(R), given an expansion his- 
tory. We show that for any barotropic equation state, 
the functional form of f(R) and hence the gravitational 
action, is not uniquely determined by the expansion his- 
tory of the universe. Instead, for a given fluid one can 
always construct a class of gravitational actions that will 
have the same cosmological evolution as the EH action. 
As an example, we demonstrate this explicitly for cosmo- 
logically relevant solutions. 

BASIC FORMALISM OF f(R) GRAVITY 

The action for f(R) gravity is (with 8ttG = 1) (see e.g. 
0) 



S = J d i x^—gi y f{R)+C m ), 



(1) 



form which it follows in the standard metric approach (as 
opposed to the so-called Palatini approach) 



G>„ — Rf_ lv ^RQiiv — T^ v + T^ v , (2) 

" = T(R){\ 9 ^ f{R) ~ Rr{R) ) + 

+ f\R)^(g a »9pv-9»v9^)}). (3) 



where 



T?. 



The standard minimally coupled stress-energy tensor, 
TJH, from the matter Lagrangian, C rn in the action, Eq. 
Q, is related to T™ by = f™//'(i?). The equa- 
tions of motion arising from the action in a Friedmann- 
Robertson- Walker universe are 



and 



H 2 



2H + 3H 2 



f'(R) 



= - p 



(4) 



(5) 



where the energy density and pressure of the curvature 
fluid arc 



Pc 
Pc 



j±-{l(f(R) Rf(R)) 3HRf"(R)} 
jl^{R 2 f'"(R) + 2HRf"(R) + Rf"(R) 

~ 2 



~[f(R)-Rf'(R))}. 



(6) 



In addition to these, we also have a constraint equation 
for the curvature scalar 



R + 6[H + 2H Z 



= 



and we can write the continuity equation of the total 
fluid from Eqs. (QJ and J5J as p + 3H(p + p) = 0, where 
P = Pc + Pm/f'{R) and p = p c +p m /f'(R)- 

From the Einstein's equations we see that due to the 
presence of the R term in p c , we have a fourth order 
differential equation system for the scale factor of a(t). 
Given a gravitational theory, or given the form of f(R), 
one could hope to solve the arising differential equation. 
However, even for simple choices of f(R), such as f(R) = 
R+^/R, the resulting differential equation (see eg. 0) is 
complicated and one must resort to numerical methods. 
The numerical solutions are also problematic since one 
then needs information about the higher derivatives of 
the scale factor. 



UNIQUENESS OF THE GRAVITATIONAL 
ACTION 



It is clear that for f(R) = R, Equations 10} and © 
reduce to the standard Friedmann equations of Einstein 
gravity, since then p c — p c = 0. However, this does not 
guarantee that f(R) = R is the only choice that reduces 
to the standard Friedmann equations but in general stan- 
dard equations are reached if conditions 



Pc 
Pc 



Pn 



f'{R) 

Pm 

f'(R) 



Pm, 
Pm 



(8) 
(9) 



hold. Assuming that the continuity equation holds inde- 
pendently for the ordinary matter fluid, p m + 3H(p m + 
Pm) = 0, these two conditions are in fact equivalent. 

The fact that one can recover the standard Friedmann 
equations is a consequence of higher derivatives in the 
Einstein's equations. The general Friedman equations 
are higher order differential equations than the standard 
equations. Thus the standard solutions can also be so- 
lutions of the general equations without violating the 
uniqueness of the solutions of differential equations. One 
can step even a little further and show the standard Fried- 
mann equations supplemented with Eq. © are enough 
that both general Eqs. (@J and JSJ) are fulfilled. 

When the ordinary Einstein's equations hold, the cur- 
vature scalar can be written as R = —p m + 3p m , so that 
for any barotropic equation of state (expect for the ra- 
diation dominated case), p m = p m {Pm), the curvature 
scalar is a function of the matter density, R = R(p m ). 
One can then write, at least formally, p m — p m (R), e.g. 
for cold dark matter p m = —R. Therefore, the quantity 
HR appearing in the expression for p c can be written 
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as HR = -{p m - Ska 2 ){p m + p m )R'(pm). In a flat 
universe, the curvature fluid is hence a function of R 
only, p c — Pc(R) and therefore the condition, Eq. ||SJ), 
gives a second order homogeneous differential equation 
for g(R) = -R+f(R): 



S^(p m +Pm)g"(R) - (^R + Pm )g'(R) + l -g{R) = 0, 

(10) 

where p m and p m are given as functions of R only. Hence 
for any barotropic fluid, there exists a general f(R) with 
two integration constants (expect if p m + p m — which 
leads to a first order differential equation) that will lead 
to the exactly same background evolution as in the stan- 
dard EH gravity. In other words, for any barotropic equa- 
tion of state f{R) gravity constructed in this way, will 
always have among its solutions the same solution as the 
solution arising from EH gravity. 

In a non-flat universe, the conclusions are unchanged 
but the technical details are somewhat more complicated 
since now one cannot write HR as a function of R only 
due to the presence of the fc-term. Instead, one needs to 
express everything in terms of the scale factor a, resulting 
in a differential equation for F(a) = f(R(a)) to be solved. 
After inverting R = R(a) one obtains an expression for 

/OR)- 



CONSTRUCTING EQUIVALENT 
GRAVITATIONAL ACTIONS 

To illustrate the general conclusions of the previous 
section, we construct explicit gravitational actions that 
have identical background expansion as the solution aris- 
ing from the EH action among their solutions. 

Consider first a fluid with a constant equation of state, 
p m = wp rn , where w is a arbitrary constant, excluding 
cases w = —1, 1/3. Now p m = /9 m ,oa~ 3 ^ 1+ ™\ and as 
explained above, expressing all in terms of R, writing 
g(R) = G(u), where u = R/(p m .o(3w — 1)), we obtain 
the general equation for G: 
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1 



3w- 1 
3w + l 
3w- 1 



1-ku 3 <™+ 1 ) 



uG'(u)+G(u), 



'G"(u) 



(11) 



where k = k/p m . In the case of a flat universe, k = 0, 
Eq. (|11|) can be solved: 



j(R) = c+i-R)^^ + c_(-R) a - 



O) 



(12) 



where a±(w) = (9w + 7 ± y/9w 2 + 78w + 73)/(12(w + 
1)) and c± are constants (note the R is negative with 
our conventions). In contrast to the k ^ case, the 
spatially flat solution does not depend on the (present 
day) density p m ,o- In the special case of a flat matter 



dominated universe, w = 0, k = 0, the solution translates 
to 



f(R) = R + c + (-Rf+^/ 12 
+ c_(-i?)( 7 -^>/ 12 . 



(13) 



Note how the other solution actually grows with \R\ i.e. 
its effect is larger at early times. 

The radiation dominated case w = 1/3 is not included 
to the solution given above because now R = and the 
scale parameter can not be expressed in terms of R. How- 
ever, this means that for radiation domination R = and 
therefore from Eqs JSJl and any function g(R) with 
g(0) = g'(Q) = will give the same radiation dominated 
expansion. The other excluded value w = —1, corre- 
sponding to the deSitter model, gives R = — 4pA- Hence 
R = and Eq. JTDJ simplifies to Rg'(R) - 2g(R) = 
giving g{R) — cR 2 . This solution also coincides with the 
limiting value of a + -solution: a + (w — > —1) — » 2. 

Another interesting example is the flat ACDM model 



a(t) 



1/3 



sinh 2 / 3 ^ 



(14) 



in which case the differential equation for g(R) reads as 

(R + 3/3) (iZ + W(R) - (\r + P)g'{R) - l -g{R) = 0, 

o o 

(15) 

where we have defined (3 = 3£IaH 2 . This equation 
has a solution in terms of the hypergeometric function 
F(a, (3; 7; z) and the solution in terms of f{R) is 

f(R) = R + c + z a +F(-a + ,l-a + ;a_-a+ + l;-) 

2 z 

3 1 

+ c_z a -F(-a_,--a_;a + -a_ + l;-) (16) 

where z = —3 — R/f3 and a± = a±(w = 0). It is easy to 
see that in the limit R — > — 00 (or z — > +00) we recover 
the matter dominated case l|13|l as expected. 



STABILITY 

An important question in considering f(R) models as 
realistic theories of gravity is the stability of the ground 
state 0|- The stability criterion can be formulated as a 
condition on the sign of the potential given by 



U(R) 



1 fW(R)fW{R)R _ f w {R)R 
3 3/( 2 ) (R) 2 3/( 2 ) (R) ' 

fM(R) 2f(R)fW(R) 
3/( 2 )(i?) 3fW(R) 2 



(17) 



where we have assumed that the universe is homogeneous 
(note that we have different metric conventions compared 



4 



to 13]). If U(R) > for classical, unperturbed solution, 
then the linear perturbations 6R are oscillatory without 
exponentially decreasing or decaying modes. In terms of 
function g(R) the stability condition obtained from Eq. 
(|T7jl reads as (for g" ^ 0) 



(2g-Rg' 



,(3) _ 



1 



g'-Rg")i 



> o. 



(18) 



Equation (|18() shows that for the radiation dominated 
case gravitational stability is achieved simply if g"(0) < 
and for the deSitter case (where g = cR 2 ) if c < 0. In 
order to have an understanding whether stable choices 
of f(R) exist in the more general case where p = wp, 
we consider the two solutions f±(R) = R + c±(— R) a± 
separately and leave more general considerations for fu- 
ture work. We find that both solutions are stable in two 
separate regions, 



c± < 
c± > 



6(1 + w)(- J R) 1 - Q±(lu) 



17+ 15w T V9w 2 



78w + 73 



(19) 



and hence one can always construct stable extensions of 
the EH action for such a fluid. 



CONCLUSIONS AND DISCUSSION 

In the present letter we have considered the uniqueness 
of the gravitational action among f(R) gravity models. 
We have showed that for any barotropic fluid, the evolu- 
tion of the scale factor of the universe does not uniquely 
determine the form of f(R), and hence nor the grav- 
itational action. Instead, one can always construct a 
f(R) theory that will have the same background evo- 
lution among its solutions. In other words, cosmology 
covered by standard Friedmann equations has always a 
generalized counterpart with the very same classical evo- 
lution independently of the actual equation of state of 
the ordinary matter. 

As an example, we have constructed explicit forms of 
f(R) in a number cosmologically interesting cases. By 
considering small fluctuations around the ground state, 
we have shown that such modified theories of gravity are 
also stable. 

In this letter we have not directly addressed the dark 
energy problem. However, our results are likely to be 
relevant to such considerations. Assuming that we have 
a matter dominated universe, we now know the form of 
f(R) we need in order to reproduce the matter domi- 
nated regime exactly. This will guide us in considering 
choices of f(R) that will lead to late time acceleration 
while preserving standard matter dominated expansion 
at early times. 

These considerations are also closely related to the 
models obtained by conformal transformations. Al- 



though the gravitational actions constructed in this ar- 
ticle may be rather unintuitive, their counterparts in 
scalar-tensor -gravity may seem natural. Moreover from 
the analysis presented here we know that some highly non 
trivial, probably non-minimally coupled scalar-tensor 
theories reproduce exactly the standard evolution of the 
background metric. 



Since the background expansion alone cannot distin- 
guish between different choices of f(R), one must study 
perturbations in such cosmologies. This will ultimately 
guide us in choosing the correct gravitational action 
which is possibly not the simple Einstein-Hilbert action. 
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